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Abstract: The existence and detection of scalar fields could provide solutions to long- 
standing puzzles about the nature of dark matter, the dark compact objects at the centre 
of most galaxies, and other phenomena. Yet, self-interacting scalar fields are very poorly 
constrained by astronomical observations, leading to great uncertainties in estimates of the 
mass and the self-interacting coupling constant A of these fields. To counter this, we 
have systematically employed available astronomical observations to develop new con- 
straints, considerably restricting this parameter space. In particular, by exploiting precise 
observations of stellar dynamics at the centre of our Galaxy and assuming that these dy- 
namics can be explained by a single boson star, we determine an upper limit for the boson 
star compactness and impose significant limits on the values of the properties of possible 
scalar fields. Requiring the scalar field particle to follow a collisional dark matter model 
further narrows these constraints. Most importantly, we find that if a scalar dark matter 
particle does exist, then it cannot account for both the dark-matter halos and the existence 
of dark compact objects in galactic nuclei 
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1. Introduction 

A large number of astronomical observations little doubt about the existence of dark matter 
(DM) and its dominant role in the matter composition in the Universe [jlp. An even larger 
number of candidates has been proposed to play the role of this unknown component. 
Among those candidates, ultra-light scalar fields have been extensively studied, providing 
a robust paradigm for DM particles [Q, [|, |5p and dark energy [0]. Another scalar field, 
the inflaton, is a crucial ingredient in inflationary cosmology [[7L H, §]. The success of 
cosmological and astrophysical models that invoke scalar fields has, in addition, motivated 
further studies of self-gravitating systems made of scalar fields, such as non-topological 
solitons [0], oscillatons 0, P, 0] and boson stars (BS) |5j. 

In particular, when considering a BS the scalar field is associated with a spin-zero bo- 
son with mass and self-interacting coupling constant A. The essentially unconstrained 
freedom in choosing the free parameters characterizing the scalar field (namely its mass 
and self-coupling constant), allows one to build BS-models that can account for many 
different astrophysical objects: from galactic dark matter halos [jT(5|, [TTj, [TS], [T9p, to dark 
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compact objects (DCOs) like MACHOS $Q [HJ] and black hole candidates [@ ||, 0]. 
Clearly, the only way to restrict the very large space of parameters spanned by BS models 
is by exploiting astronomical observations and use them to set constraints on the properties 
of the scalar field. This is, in essence, the goal of this paper. 

The recent advances in high-angular resolution instrumentation have provided the pos- 
sibility to study the central regions of galaxies with unprecedented precision. Space-borne 
telescopes, such as the Hubble Space Telescope, and ground observations using adaptive 
optics have opened a new window on the innermost central regions of galactic nuclei. Par- 
ticularly striking is the present ability to study the kinematics of stars or gas in regions of 
milliparsec scale for the Milky Way [|25|, ^(J and of sub-parsec scale for external galax- 
ies [ p7| , ^8], |29|]. One of the most intriguing conclusions of these observations is that dark 
compact objects (DCOs), with masses ranging between M DCO — 10 6 — 10 9 M , are hosted 
at the centre of most non-active galaxies. 

In the case of our own galactic centre, the study of the innermost stellar dynamics has 
provided very convincing evidence for the existence of a DCO associated with the radio 
point source Sagittarius A* (Sgr A*). The closest Keplerian orbits examined are those of 
the so-called "S-stars" (also referred to as "SO-stars"). By following the orbit of one of 
them, i.e., the star S2 (also referred to as S02), the mass of Sgr A* has been estimated to 
be about 3.7 x 10 6 M Q within a volume with radius no larger than 6.25 light-hours [|2~6[ |3Q|]. 
More recent data based on 16 years of observations has reduced the uncertainty and set the 
mass of the central DCO to ~ 4.1 ± 0.6 x 10 6 M [|1],]32|, Hi- 
Assuming that the DCO in Sgr A* is a BS, we have used the observational restrictions 
on the volume and mass to set substantial constraints on the mass of the scalar field and on 
its self-interacting coupling constant. The approach followed using the observational data 
for Sgr A* can be applied to any other galaxy for which there is an accurate measurement 
of the mass and dimensions of the central object. As a result, we have also considered 
the case of a nearby galaxy, NGC 4258, obtaining additional constraints, which set much 
tighter limits when considered in combination with those of Sgr A*. 

Finally, we have considered several DM candidates and discussed how present astro- 
nomical observations can be exploited to set additional limits on the scalar-field properties 
invoked by these DM models [Q, [16], [17|, [18|, [jj^]. While some of these constraints over- 
lap with those obtained when considering BSs, others cover a distinct region of the space 
of parameters and thus imply that a scalar field that could explain the rotation curves in 
nearby galaxies cannot be the same composing a BS representing the DCOs in galactic 
nuclei. 

The structure of the paper is the following one. In Sect. ^ we briefly review the prop- 
erties of BSs in the two regimes of weak and strong self-interaction. Section |3| is instead 
dedicated to the issue of the maximum compactness of a BS, where we show that the def- 
inition of an effective maximum compactness for a BS is not only possible, but actually 
useful when associating BSs to DCOs in galactic centres. In Sect. |] we therefore derive 
our constraints on the scalar-field properties using both considerations on BSs models and 
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of DM models. Finally, our conclusions are contained in Sect. |5|. 



2. Boson star generalities 

Boson stars are solutions to the Einstein Klein-Gordon system of equations 

G, u = torGT^ , (n-^j=0, (2.1) 

where □ = (1/ y/—g)d ll [s/^gg IM 'd v ], We have adopted natural units h — c — 1, and hence 
G = l/m 2 p , where m p is the Planck mass. The stress energy tensor T^ u corresponds to a 
complex scalar field minimally coupled to gravity and with a scalar potential V(|$| 2 ), 

Tfi V = \[d^*d v (Sf + M $0„$*] - i[$ *' Q $ Q + Vm 2 )] • (2-2) 

We will restrict ourselves to the case where the potential of the scalar field is given by V = 
|m||$| 2 +||$| 4 , where is the mass of the scalar field and A its self-interaction coupling 



constant. Originally BSs were introduced by Kaup [ |34| ] and studied later by Ruffini and 



Bonazzola [p5Q, who considered just the mass term in the scalar potential. The case of 



self-interacting scalar fields were later discussed by Colpi and collaborators in [p6|]. If 
for simplicity the scalar field is assumed to be spherically symmetric and with a harmonic 
time dependence, $(r, t) = <j)(r)e~ luJt , the space-time part is static. Considering spherical 
symmetry, the line element can be written as ds 2 = —B{r)dt 2 + A(r)dr 2 + r 2 dVL 2 and 
the Einstein Klein-Gordon equations reduce to a system of ordinary differential equations 
for the metric functions A and B = rr&B / 'co 2 , and for the scalar field <fi [^6|]. Additional 



quantities appearing in the equations are o = y?Aix /m? p (p, which reduces to the density in 
the Newtonian limit, and the dimensionless self-interacting coupling constant 



Am 2 

A = - — % • (2.3) 

The freedom of choice of A and is reflected in the freedom of choosing A and 
it has been shown that even when A < 1, the resulting configuration may differ signif- 



icantly from the non-interacting case [[36|]. For sufficiently small A, on the other hand, 
self-interaction may only be ignored if A m|/m 2 , that is, if A 1. On the other hand, 
configurations with A ^> 1 are interesting as their masses may be comparable to those of 
their fermion counterparts if A ~ 1. This illustrates the relevance of taking both regimes, 
A ~ 1 and A 3> 1 in the analysis of BS configurations. In the next two subsections we 
review them briefly for completeness and to introduce a notation that will be useful later 
on. 
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Figure 1: Gravitational mass M99 of equilibrium configurations of BSs for different values of a c 
and A in the regime A3 <C 1; black squares indicate configurations with the maximum mass. The 
left panel, in particular, shows a three-dimensional view of the space of solutions, while the right 
panel offers a view of three slices given values of A. 

2.1 Weak self-interaction (A 3 < 1) 

We next we derive equilibrium configurations for the Einstein Klein-Gordon system in the 
weak self-interaction regime, i.e., when A 3 = A/1000 <C 1, following a procedure similar 
to the one discussed before in [|T4|, O, 123,123]]. In particular, in this case the Einstein Klein- 



Gordon system of equations can be solved as an eigenvalue problem for the radial metric 
function at the BS' centre B c after specifying a finite but arbitrary value for the "central 
density" o c . Additional boundary conditions are those of regularity at the centre, i.e., A c = 
1, a' c = 0, where the prime denotes a derivative with respect to the new dimensionless 
radial coordinate x = rm^, and of asymptotic flatness, i.e., o(x = 00) = 0. 

The system is solved numerically for different values of a c and A using a standard 
shooting method (given a a c , there is a unique value of B c for which the boundary con- 
ditions are satisfied). Because we use a finite numerical domain, the condition for a at 
infinity is in fact demanded at the outermost point of the domain x out , and the shooting 
procedure is performed for different values of x out . As x out increases, the shooting param- 
eter B c converges, and we choose the solution by requiring that the condition o-(x out ) = 
holds within a prescribed tolerance (i.e., o~(a: out ) < 10- 9 ). 

Once the equilibrium solution is obtained, its mass is computed simply as 

M = ^ ( 1 - — i— ^ . (2.4) 



^out) 

We note that since the scalar field a decays exponentially, the radius of the star is, at least 
formally, not finite. In practice, however, we define Afg 9 = 0.99 M, that is as 99% of the 
total gravitational mass M; as a result, we also define R 99 as the radius containing M 99 
and set this to be the "effective radius" of the BS. As a final remark we note that in view 
of the units adopted, both quantities M and R are dimensionless variables. The physical 
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Figure 2: Gravitational mass M99 of equilibrium configurations of BSs for different values of a c 
in the regime A3 > 1. Also in this case the black square indicates the configurations with the 
maximum mass. 

units for the mass M and radius R can be recovered by using the following relations 

m 2 * h 

M = M — - , R = R — . (2.5) 

Figure [I] reports the gravitational mass M 99 of equilibrium configurations of BSs for 
different values of a c and A in the regime A 3 <C 1. In particular, the left panel shows a 
three-dimensional view of the space of solutions, while the right panel offers a view of 
three slices given values of A. In both panels, black squares indicate configurations with 
the maximum mass, thus distinguishing stable solutions (to the left of the squares) from 



unstable ones (see discussion in Sect. 2.3). 



2.2 Strong self-interaction (A3 > 1) 

The solution of the Einstein Klein-Gordon equations is much simpler when considering 
instead the strong self-interaction regime, i.e., when A 3 3> 1. In this case, in fact, intro- 



ducing the new dimensionless variables [36] 



<r„ = A 1/2 (T , (2.6) 
Xif = A~ 1/2 x . (2.7) 

and neglecting terms of order 0(A _1 ), the Klein-Gordon equation can be solved alge- 
braically to yield p^[ ] 

I \ V2 



a * = \B~ 1 ) ' (2 ' 8) 

so that B c = l/(a 2 c + 1). Furthermore, the full set of the Einstein Klein-Gordon equa- 
tions does not depend on A at first order and can be solved much more easily. As in the 
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a a. 

Figure 3: Effective compactness C = Mgg/Rgg for a spherical BS shown as a function of the 
central density a c . Shown in the left panel are the values in the regime of weak self-interaction 
A3 <C 1, while the right panel reports the values in the regime of strong self-interaction A3 3> 1. 
Also in this case the black squares in both panels indicates the configurations with the maximum 
mass. 



weak self-interaction regime, equilibrium configurations are obtained solving a set of ordi- 
nary differential equations with the same boundary conditions. One important difference, 
however, is that in this limit the density a does go to zero at a finite radius x s and, as 
for ordinary fluid stars, the gravitational mass M is defined by Eq. ( [2.40 with x ou t = ^s- 
In analogy with the previous treatment, we define the effective radius of the BS as Rgg, 
namely as the one containing 99% of the total gravitational mass 1 . Finally, because the 
introduction of the new variable 2*, see eq ( ]2!6|) , the physical units can be recovered by 
using the following relations 

M = A 1/2 M^ , R = A 1/2 R — . (2.9) 

Shown in Fig. |2| is the gravitational mass M 99 of equilibrium configurations of BSs for 
different values of a c in the regime A 3 3> 1 . Also in this case the black square indicates the 
configurations with the maximum mass, thus distinguishing stable solutions (to the right 
of the square) from unstable ones (see discussion in Sect. Of). 



As a final remark, we note that a general analytic expression for the mass and com- 
pactness of a BS in the strong self-interaction regime can be given and that fits well the 
data presented in the right panel of Fig. |2|. Such expressions are useful as they allow one 
to obtain simple estimates to be used in astronomical observations. 

2.3 Stability of Boson Stars 

Clearly, the stability of BS configurations is a basic requirement when considering them as 
suitable models of astrophysical objects. Numerical [ |37| , [38| |39| ] and analytical p0| , |4T| , [421 ] 



'Note that since in the case A 3 3> 1 the BS has a finite radius, the definition of i? 9 g is sttictly speaking 
not necessary, but it is nevertheless useful to maintain consistency with the case A 3 <C 1. 
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stability studies of BSs agree on a general result: nodeless BS configurations are stable 
under finite radial perturbations if their scalar field central value a c is smaller than a critical 
value cr cr i t ; configurations with a c > cr crit , in fact, are either collapse to a BH or disperse 
at infinity @. 

Note that, for any given value of A, the critical density cr crit also marks the maximum 
mass corresponds the maximum mass of the BS, i.e., 



M max (A) = M{a c 



(T, 



crit) 



A). 



(2.10) 



In the regime of A 3 <C 1, the maximum mass clearly increases as A increases as 
shown in the two panels of Fig. [I]. On the other hand, when A 3 1, one can exploit the 
fact that the system of the Einstein Klein-Gordon equations can be identified with that of 
a perfect-fluid star for which the effective equation of state is 



V 



1 + 



3p 
4p 



1/2 



- 1 



where 



(2.11) 



(2.12) 



and p = mi/ A. Note that in the limit p <^ p (which is equivalent to P <C p i.e., the 
Newtonian limit) Eq.( |2.1 1| ) reduces to P = p 2 / (16p )> thus representing the equation of 
state of an n = 1 polytrope. As a result, the well-known theorems of stability of fluid stars 
may be applied to determine the stability of the BSs. In particular, BSs with 1/B C smaller 
than a critical value 1/B crit are stable, while those for which 1/B C > 1/B crit will either 
collapse to a black hole or disperse to infinity. In practice, when A 3 ^> 1, the maximum 
mass is M max ~ 0.22 and is attained at <T cri t ~ 0.97; this is shown with a black square in 
Fig- 1- 



3. Maximum Compactness for a Boson Star 

Since we are interested in modelling astrophysical DCOs as stable BSs, it is important to 
determine whether or not a BS admits a maximum compactness, just as it is the case for an 
ordinary relativistic star. At first it may appear that this question is not even well posed as 
BSs are in principle infinitely extended and hence always with a zero compactness. How- 
ever, even an infinitely extended BS is in practice centrally condensed, i.e., with a finite 
effective radius _R 99 . Hence, it is mathematically legitimate and astrophysically reasonable 
to consider whether the "effective compactness" of a BS, defined as 

C{a c ,A) = , (3.1) 

-K99 
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and thus measuring the compactness of a BS with mass M 99 in a radius Rgg, can be nonzero 
and is even upper bounded. As we will show in what follows, this is indeed the case and 
C(a c , A) < 0.16 for all values of A. 

Before entering the discussion of the different regimes of self-interaction it is interest- 
ing to note that if computed in terms of the physical values for the effective mass Af 99 and 
effective radius i? 99 , the compactness C is independent of [(cf. Eqs ( |2.5[ ) and (|2.9|)]. 

In the case of weak self-interaction A 3 1, the compactness ( |3.1| ) is shown in the left 
panel of Fig. |]for different values of A. Note that the compactness is a growing function 
of <t c , but also that it has a local maximum. The latter, however, is reached for unstable 
configurations, since it refers to BSs whose central density cr c > a CT i t . Shown in fact with 
black squares are the configurations with the maximum mass and these are all to the left of 
the local maxima: these models mark therefore the BSs with the maximum compactness 
for a given choice of central density and coupling constant, i.e., C max (a c , A). 

A very similar behaviour is shown by the compactness in the right panel of Fig. [3J, 
which however refers to the regime of strong self-interaction A 3 ^> 1. The obvious dif- 
ference in this case is that the compactness does not depend on A (and hence on the free 
parameters of the scalar field and A) and attains a maximum value for stable configu- 
rations Cmax A^ = 0.158496114 at a crit ~ 0.97; this is shown with a black square in the 
right panel Fig. |] 2 . Interestingly, this maximum compactness is close to that of a typical 
neutron star with mass M N $ = 1-4 M & and a radius of -Rns = 12 km. 

Comparing the two panels of Fig. it is easy realize that the maximum compactness 
increases with A and has as asymptotic value for A — > oo the one given by C^^x^ . This 
is clearly shown in Fig. 01 where we report the maximum compactness C max as a function 
of the of coupling constant A. Shown with a dashed line is the asymptotic maximum 
compactness for a spherical BS (7^^^ ~ 0.158, while shown with a dotted and dot- 
dashed lines are the maximum compactnesses for a spherical star C maXiStar = 4/9 [Ql and 
the compactness of a black hole C = 1/2. 

Overall, Fig. |] shows two important results. The first one is that although BSs are 
in principle infinitely extended they can be considered to have an effective compactness 
which is upper limited. The second one is that such a maximum compactness is smaller 
than the corresponding one for a fluid star and, of course, of a black hole. To the be best 
of our knowledge neither of these two results was discussed before in the literature. 



4. Constraining and A 

In what follows we discuss how to make use of the astronomical observations of galactic 
centres (Sgr A* and NGC 4258) and of present DM models to constrain the space of pos- 
sible parameters for the scalar field mass and self-interaction coupling constant A. It 

2 Note that although C max ^ x does not consider the whole mass, the corresponding compactness 
0.158496114 is very close to the maximum one obtained when considering the whole mass and the whole 
radius and which is given by M^/Rs = 0.159807753. 
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Figure 4: Maximum compactness C max shown as a function of the of coupling constant 
A. Indicated with dashed line is the asymptotic maximum compactness for a spherical BS 
Cmax.Aoo — 0.318, while a dotted and dot-dashed lines are the maximum compactnesses for a 
spherical start C max , s tar = 4/9 and the compactness of a black hole Cbh = 1/2, respectively. 

is however important to recall that present constraints on rrij, and A are extremely loose. 



Previous works that considered BSs as DCOs candidates Q23| , [241 ], in fact, neither or A 
were related to an existing particle or particle candidate. In Ref. [[22|], on the other hand, 
some first estimates were presented but not using the astronomical observations considered 
here. 

4.1 Constraints from Sgr A* 

A first set of constraints can be set by exploiting the observations about the dynamics of 
the S-stars around Sgr A* and in particular of the smallest measured periastron for this 
group of stars. This is attained for the star S2, whose periastron is R S2 — 17 light-hour 
~ 5.9 x 1CT 4 pc. If a DCO is indeed hosted at the centre of Sgr A*, then its radius cannot 
be larger than i?g 2 or it would perturb considerably these orbits. Hence its minimum 
compactness is 

R S2 "3015' 1 } 

where the second equality has been obtained using the measured mass of Sgr A*, which 
has been estimated to be M Sgr A * = 4.1 ± 0.6 x 10 6 M Q P5[] 3 . Note that when expressed 



in terms of the Schwarzschild radius corresponding to Sgr A*, i.e., Rs = 2M Sgr A«, the 
periastron of the S2 star is R S2 — 1520 R s . 



3 We note that larger compactnesses (indeed as large as the ones corresponding to a black hole) have been 
suggested for C m i n on the basis of radio observations at wavelengths of 3.5 mm and 7 mm |47| ]. The 
measurements are particularly complex and may be contaminated by instrumental errors which are difficult 
to remove completely. In view of these uncertainties we prefer to use the more conservative but also more 



accurate estimate (4.1 ). 
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Figure 5: Constraints on the possible values of the mass and coupling constant of the complex 
scalar field. The left panel refers to the weak self-interaction regime of A3 <C 1, where the dashed 



line refers to the condition (4.3), while the solid one to the condition (4.4). The shaded region 



is clearly the union of the inequalities. The right panel is the same as the left one, but refers to 



the strong self-interaction regime of A3 <C 1, so that where the dashed line refers to (4.5), while 



the solid one to ( |4.6| ). In both panels shown with the insets are the same data but in terms of the 



dimensionless self-interacting coupling constant A. 

The observational evidence that Sgr A* is a black hole is indeed very convincing, but 
we here assume, for the sake of argument, that it is rather a BS of mass M SgrA *- In this 
case, its compactness could be in the range 

3.32 x 1(T 4 ~ C min < C BS < C max ~ 0.158 . (4.2) 

In other words, by changing and A it is possible to construct infinite BS models that 
would be compatible with the observations of the S-stars and with radii between R S2 and 



R mm = M Sgr A * /C max . In this way, the condition ( pk2| ) can be used effectively to constrain 
the space of parameters of potential scalar fields. 



The procedure adopted in practice to enforce the condition QOD is somewhat involved 
and, as in the previous Sections, we will consider separately the cases of weak and strong 
self-interaction. In particular, when A 3 1 the maximum value of compatible with the 
mass of Sgr A* can be set by recalling that the relation between physical and non-physical 
masses ( ]23| ) states that the scalar-field mass is proportional to the BS mass for any given 
value of the physical mass. Hence, after fixing M SgrA *, the mass of the scalar field will 
have to satisfy the inequality 

, M max (A) 2 
m<f> < -rr m p . (4.3) 

JWSgr A* 

The condition ( |Q1 ) when the equality holds is shown as a dashed line in the left panel of 

Similarly, a lower limit on can be found by determining, for each A, the central 
density a* and the radius R(cr*) of the BS model having the minimum compactness, C min 
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(we recall that the compactness is a monotonic function of C). Defining then the minimum 
mass as M min (A) = C mm -R(<T*), we can constrain the mass of the scalar field to satisfy the 
inequality 



. M min (A) 2 
> — m . (4.4) 

Ivl Sei A* 



The condition Q4-.4D when the equality holds is shown as a solid line in the left panel of 



Fig. |5|, while the shaded region is the union of the inequalities ( |4.3D and (|4-.4|). 



In the regime of strong self-interaction A 3 ^> 1, on the other hand, we employ a 
technique similar to the one discussed above to obtain again upper and lower limits, with 
the difference that in this regime there is only a single mass curve (cf. Fig. ||). In this case, 
the maximum is 



/O 22 m 3 

m < m™ ax = x I - £ A 1/4 ~ 2.9 x 10 5 A 1/4 eV . (4.5) 

M SgT A* 



The condition ( J4.5[ ) when the equality holds is shown as a dashed line in the right panel of 
Fig. |5|. In addition, also when A 3 3> 1, C min allows us to derive the lower limit of with 
the difference that the variable parameter now is just B c and we do not need to iterate also 
on A. As a result we obtain 

> mj n ~ 3.7 x 10 4 A 1/4 eV . (4.6) 
The condition (|4.4|) when the equality holds is shown as a solid line in the right panel of 



Fig. |5|, while the shaded region is the union of the inequalities ( |4.3| ) and fl4.6|). The overlap 



in the constraints between the two regimes is very good and hence the results obtained for 
A 3 ^> 1 can be easily extrapolated to the regime A 3 1 by using fitting expressions. 

4.2 Constraints from NGC 4258 

The procedure discussed above for the Galactic centre can be exploited in principle for 
any other galaxy for which there is an accurate measurement of the mass and dimensions 
of the central object. As an example, we here consider also NGC 4258, which is a spiral 
galaxy at a distance of about 73 — 83 Mpc. In this case, the mass of the central DCO 
has been estimated to be M N gc4258 = 38.1 ± 0.01 x 10 6 M , while the observations of 
the rotation curve require a central density of at least 4 x 10 9 M pc -3 , which implies a 
maximum radius R max ~ 36000 R s [|48|]. 

Restricting our attention to the regime A 3 1, we show in Fig. ^| as a forward- shaded 
area the upper and lower limits for (m^, A) as derived from the observations coming from 
NGC 4258. Similarly, in the regime A 3 ^> 1, the upper limit for the scalar-field mass is 
obtained by setting the maximum mass < m™ ax = ^0.22 m 3 p /M NGC 425sA 1//4 , while 
the lower limit comes from the minimum compactness C min = 1/72000. The combined 
constraint can then be expressed as (not shown in Fig. |^) 

6 - 3eV ^^^ 9 - 6xl ° 4eV - ( 4 - 7 ) 
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Figure 6: Constraints for and A derived from the observations of NGC 4258 in the limit A3 <C 
1 (forward-shaded area). Shown as a double-shaded area is the combination of the constraints 
coming from Sgr A* and NGC 4258. Shown with the inset is the same data but in terms of the 
dimensionless self-interacting coupling constant A. 



Combining these constraints with those derived in the previous Section for Sgr A*, we 
obtain the following global constraints on the parameter space obtained when modelling 
the DCOs as BSs 



3.7 x 10 4 eV < 



m0 < 9.6 x 10 4 eV . 



A 1 / 4 



(4.8) 



The corresponding region is shown as a double shaded region in Fig. ^| and clearly sets a 
tighter constraint on the possible values of and A. 

4.3 Constraints from Dark-Matter Models 

As mentioned in the Introduction, constraints on the properties of the scalar field can be 
imposed also exploiting cosmological observations. We recall that the Lambda Cold Dark 
Matter model (or ACDM as it is usually referred to 4 ) is a standard model of big bang 
cosmology, which attempts to explain within a single framework a number of cosmolog- 
ical constraints and observations. These are, for instance, the existence and properties of 
the cosmic microwave background, the abundances of light elements (hydrogen, helium, 
lithium, the large-scale structure of the universe in terms of galaxy clusters, and, more 
recently, also the accelerating expansion of the universe observed in the light from distant 
galaxies and supernovae. 

In spite of its many successes, some of the predictions of the ACDM scenario are in 



contrast with observations. One example is the so-called the "core/cusp problem" [|4J| p0| ], 
where numerical N-body simulations assuming a ACDM cosmology predict the generation 
of diverging cusps in the central density of the DM distribution in galaxies, and which are 
incompatible with observations at galactic scales [|5T|]. Another example is the so-called 



4 Note that the A here is the cosmological constant and should not be confused with the self-coupling 



constant defined in (2.3) 
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the "dwarf-galaxy problem" [|52], where again the ACDM predicts an overproduction of 
dwarf galaxies in the local group, which is not corroborated by observations. Some pos- 
sible solutions to these problems have been suggested even within the ACDM paradigm. 
For instance, the core/cusp problem could be alleviated if dynamical friction [ |53"| ] or stellar 
feedback [ |53| ] are included. Similarly, for the dwarf-galaxy problem, it was suggested that 
accretion of gas onto low-mass halos is not enough to make them observable. Possible 
reasons of this inefficient accretion could be supernova feedback [155], |5H|, photoioniza- 



tion D57| ] or reionization [|58fl. It has also been suggested that this may simply be a problem 
with our ability to detect the missing satellite galaxies in the Local Group, or that these are 
in a region of the sky that has not been surveyed yet []59|]. Alternative solutions suggest 
instead that DM might be treated as extremely light bosonic (dark) matter. Some models 
based on this idea are the so-called "fuzzy dark matter" [Q], or the scalar- field (dark) mat- 
ter [JJ, |], 0| . Both suggestions predict flat density profiles at the centre of galaxies [ |60| ] 
and fit the abundance of dwarf galaxies [0, 0, §]. An alternative suggestion is that the 
DM particle could have a very strong self-interaction, but negligible annihilation [^TJ . In 
this case, the core/cusp problem is alleviated since the strong self-interaction increases the 
scattering among DM particles. The scattering cross section is then so high that collisions 
at the centre of galaxies are very so frequent that dark matter particles scatter out the centre 
as fast as they are accreted, thus effectively preventing a growth in the density. 

In what follows we briefly summarize the main results of the scalar-field DM and 
the "collisional" DM scenarios and discuss how to use the present observations to further 
constraint the space of parameters for a complex scalar field. We recall that the collisional 
DM model of [ |5T] ] is also referred to as "self-interacting" DM because the scattering cross 
section is much larger than that of standard WIMPS, that collisions among DM particles 
are very frequent. On the other hand, the prototype of a DM candidate consisting of ultra- 
light scalar field is the axion, which has been introduced as a solution to the strong CP 
problem [ |62| , |63| ]. The mass of the axion has been estimated to be in the range 1CT 5 eV < 
m a < 10~ 3 eV []64t| . However, in order to explain the singular galactic core and suppress 
the formation of dwarf galaxies, the axion's mass should be orders of magnitude smaller. 
In the absence of self-interaction, the mass should in fact be ~ 1CT 22 — 10 _24 eV, 
which then gives a reasonable fit to the rotation curves [(H], |T7| ]. 

At the same time, the cosmological study of the scalar-field DM performed in [0, |55J, 
Wj \ for different potentials V(<f>) have shown to be able to reproduce all the features of the 



standard ACDM in the linear regime of perturbations. Furthermore, a DM model with a 
complex scalar field having a quartic potential was studied in [fT^] , and the fit of rotational 
curves for dwarf galaxies was shown to be robust provided that 

2-7eV<^<2.9eV. (4.9) 

This is shown as a thick solid line in Fig. [7]. 

A different suggestion has been made in Ref. [^TJ], where the DM is proposed to be 
cold, non-dissipative and self-interacting with a scattering cross-section per particle mass 
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Figure 7: The constraints on and A coming from the different DM models. Shown with a 
dashed magenta area the allowed values for coming from the collisional DM model, while shown 
with a black region are the corresponding values as constrained from the measurement of the DM 
halo. 

given by 

^ = 10- 25 - 1(T 23 — , (4.10) 

GeV 

and independent of the particle itself. This model solves the small scale problems of 
ACDM and N— body numerical simulation confirms this W7ft. Since the conjecture is 
made without assuming a particular model for DM, we can derive the values for (m^, A) 
for a scalar field with self-interaction. For that, we introduce as the interacting potential 
V = A0 4 /4. Thus, the scattering matrix element is M.{4>4> 00) = ^A and the resulting 
cross-section is 

A 2 A 2 

(7(00 -> 00) = (T H = — = — J > C 4 - 11 ) 

since the square of the centre of mass energy is s = 4m 2 . By requiring that the cross- 
section of Eq. ( [4.1 1[ ) has the strength of the collisional DM of Eq. fl4.10| ), we find that 



must be 

9.5 x 10 5 eV < — £ < 9.5 x 10 7 eV. (4.12) 

The constraints coming from the different DM models are summarized in Fig. [7], where we 
show with a dashed magenta area the allowed values for and A coming from the colli- 
sional DM model [cf. eq. ( |4. 10[ )], while shown with a black region are the corresponding 



values as constrained from the measurement of the DM halo [cf. eq. (|4.9|)]. 
4.4 Collecting all constraints 

We can combine all of the constraints derived so far and conclude that a scalar field which 
can account for DM candidate and fulfill the constraints imposed by the observations on 
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Figure 8: Constraints on (m^, A) based on available stellar kinematic data for Sgr A* (narrow 
upper stripe filled with solid blue lines), as well as for NGC 4258 (broad lower stripe filled with 
dashed orange lines) and their intersection (middle stripe filled with black crosses). The solid 
magenta stripe that crosses all other stripes is derived from the requirement of collisional DM and 
demarcates the combined restrictions (red solid region). We also show the limits on the parameters 
derived from the models of galactic DM halos based on scalar fields (solid black line). 

the DCO at the centre of our Galaxy and NGC 4258 must have a mass and a self-interacting 
coupling constant in the ranges 

6.5 x lCT 9 < A < 4.2 x lCT 3 , 

332 eV < m < 2.46 x 10 4 eV. (4.13) 

This region is marked as "combined" in Fig. || which reports also the constraints coming 
from modelling the DCOs as BSs and discussed in the previous Sections. 

We also note that although the lower limit of the region constrained from the observa- 
tions of NGC 4258 is not too far from the limit set by condition Q4.9|), the latter is at least 
four orders of magnitude smaller than the region representing the overlap between colli- 
sional DM models and the kinematic constraints from galactic observations. As a result, 
we conclude that a scalar field that could explain the rotation curves in nearby galaxies 
cannot be the same composing a BS representing the DCOs in galactic nuclei. 
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5. Conclusions 



We have used the most accurate astronomical observations on the kinematics of galaxy 
centres to set constraints on the mass and self-interaction constant A of complex scalar 
fields, which have been employed in the literature to model BSs as DCOs in galactic 
centres or as DM candidates on cosmological scales. More specifically, we have used the 
estimates of the mass and volume for the nucleus of our Galaxy and of NGC 4258 to limit 
the possible models of BSs when these are advocated to account for such observations. 

These constraints, together with those obtained by considering possible DM candi- 
dates have allowed us to restrict the possible values of and A to the ranges 

6.5 x 1(T 9 < A < 4.2 x 1(T 3 , 
332 eV < < 2.46 x 10 4 eV , 

that is to a region that spans six orders of magnitude for A and one order of magnitude for 
nifj,. To the best of our knowledge, these constraints for the scalar field are by far tighter 
than any other discussed so far in astrophysical or cosmological context. In addition, since 
the space of parameters has two non-overlapping constrained regions, we conclude that a 
scalar field that could explain the rotation curves in nearby galaxies cannot be the same 
composing a BS representing the DCOs in galactic nuclei. 

Our analysis here has been limited to the simplest possible scenario in which the DCO 
is not a black hole but a BS. However, as long as the presence of a BS at the galactic centre 
cannot be ruled out, a more plausible configuration is one in which the BS actually hosts 
at its centre a massive black hole of comparable mass, i.e., a hybrid BS. Mixed models 
for BS have been proposed in the past, for instance when studying the boundary between 



stable and unstable equilibrium configurations of cold boson-fermion stars (see Q68| , |14| ] 
and references therein). 

While we are not suggesting that a hybrid BS hosting a BH at its centre is the most 
natural configuration, we do believe it is more appealing than a simple BS and for a num- 
ber of reasons. First, as long as the size of the black- hole horizon is much smaller than 
the scalar-field Compton length, the accretion of the scalar field onto the black hole is 
very small J69| ] and hence a stationary model can be constructed (if this was not the case 
the scalar field would eventually be all accreted by the black hole, as recently suggested 
in [J7DQ). Second, a hybrid BS automatically satisfies all the constraints usually set when 
considering the DCO as a massive black hole. Finally, by having a black hole at its cen- 
tre, a hybrid BS does not need to account for the electromagnetic emission that would be 
otherwise produced by the matter condensing and heating-up at the centre of an ordinary 
BS. 

Together with these advantages, however, a hybrid-BS model has the important draw- 
back that it may be very difficult to distinguish it from a pure black-hole solution, the 
differences being so small that they may be below the present (and possibly future) ob- 
servational limits from electromagnetic radiation. In this respect, the gravitational- wave 
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observations that will be made possible by the space-borne LISA mission, may well be 
determinant. LISA will in fact detect a large number of extreme-mass ratio inspirals or 
EMRIs (see, [ |7T| ] for a recent review) and because the associated gravitational waveforms 
probe a region of the spacetime very close to the black hole, they may work as telltale 
about the presence of a scalar field surrounding the black hole. Some work on the ef- 
fects of matter fields on the gravitational-wave emission from EMRIs has already been 
done [[72], [73p, and extending this question also to scalar fields will be the focus of our 
future research. 
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